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Abstract 

We present numerical simulations of a theory of resonant processes in a 
frozen gas of excited atoms interacting via dipole-dipole potentials that vary 
as r~^, where r is the interatomic separation. The simulations calculate time- 
dependent averages of transition amplitudes and transition probabilities for 
a single atom in a given state interacting resonantly with a uniformly dis- 
tributed random gas of atoms in a different state. The averages are over 
spatial configurations of the gas atoms, which are held fixed while the reso- 
nant interaction creates a Frenkel exciton that can travel from atom to atom. 
We check that the simulations reproduce previously known exact results when 
the exciton is not allowed to propagate [Phys. Rev. A 59, 4358 (1999)]. Fur- 
ther, we develop an approximation for the average transition amplitude that 
compares well with the numerical results for a wide range of values of the 
system parameters. 
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I. INTRODUCTION 



Frozen gases are a new and in some ways ideal laboratory to test and refine our under- 
standing of quantum theory in a complex system. With laser technology, the translational 
temperature of the gas can be lowered to the point where it can be ignored when discussing 
electronic processes, and at the same time one can manipulate and detect these processes 
in the gas with an extraordinary selectivity and precision. Frozen Rydberg gases have the 
added advantage that their states are well understood and electronic processes occur on a 
microsecond time scale, easily allowing for time- resolved spectroscopy. These gases provide 
an excellent means of studying many-body effects, because the dipole-dipole interactions 
which are important for binary collisions at room temperature are still dominant in the 
frozen gas, but now involve many atoms that are, to a first approximation, standing still. 

Pioneering experiments on resonant processes in frozen Rydberg gases have been carried 
out by Anderson et al. and Mourachko et al. ^ . The work presented here was motivated, 
in part, by the desire to understand those experiments and the subsequent work of Lowell 
et al. highlighting the dynamic aspects of these resonant, many particle systems. The 
measurements of Refs. are of mixtures of ^^Rh atoms initially prepared in the 23s and 

33s states, henceforth to be called the s and s' states, respectively. There are initially N 
atoms in the s state and N' atoms in the s' state. The transition ss' — > pp' is monitored, 
where p refers to a 24pi/2 state and p' refers to a 34p3/2 state, at and near resonance. 
What is measured is Np//N', where Np/ is the number of atoms in the p' state after an 
elapsed time t. (Measuring Np/N gives nothing really new, since Np = Np' at all times.) 
Some of the experimental features we wish to understand are the rapid rise followed by slow 
approach to saturation of the signal and the width of the resonance lineshape (in the detuning 
A = e'p + ep — e'^ — e^) as a function of time. The problem is not trivial because an initial 
ss' — > pp' transition is followed by the propagation of sp —>■ ps and s'p' p's' fluctuations, 
which are responsible for the broadening of the resonance. Although we address specifically 
the Rb experiments, the theory we present here is broadly applicable to analogous processes 
in any frozen gas of excited atoms. 

In an earlier paper |^ we developed the general formalism, obtained exact results for the 
case when resonance broadening is negligible, and introduced simple approximations for the 
broadening. The same formalism and notation are used in this paper. However, in Ref. |^ 
we introduced the down and up states of an effective spin to describe the states s and p, 
and similarly for s' and p', but in this paper we prefer to describe an s — > p transition as 
the creation of a Frenkel exciton which then propagates through the medium. In this 

language, the initial creation of a.pp' pair launches two excitons, of types sp and s'p', but we 
consider only the limit where the propagation of the s'p' exciton is relatively unimportant 
and can be neglected. We prefer the exciton language because the literature on excitons in 
solids and liquids generally considers the dipole-dipole forces to be of major importance, as 
they are in a frozen Rydberg gas, and this literature has been useful to us in developing 
approximate analytical treatments. On the other hand, in spin systems the dominant role is 
usually played by the short-ranged exchange interaction and dipole-dipole forces are given 
comparatively little attention. 

In this paper we present the results of numerical simulations for the case where A^' = 1 
and N is large, which corresponds to the situation where a single s' atom is initially present 
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and a single sp exciton is launched by the resonant ss' pp' transition. The calculated Npi 
can be directly compared to the experimentally determined Npi/N' in the limit A^' -C N. 
When exciton propagation is suppressed, exact results are obtained by the methods of Ref. 
and serve as a check on the numerical simulations. 

We also develop analytical approximations for the average amplitude of a resonant tran- 
sition and compare them with the results of the simulations. We find that the Cayley tree 
approximation agrees well with the simulations for a broad range of conditions, and is much 
better than the simple Lorentzian approximation introduced in Ref. 0. The details of the 
numerical simulations are briefly outlined in Section |T|. Analytical approximations for the 
transition amplitude and for the exciton spectral density are developed in Section Re- 
markably, explicit formulas are obtained in the Cayley tree approximation. Representative 
graphs of the results, both numerical and analytical, are shown and discussed in Section 0, 
which is followed by the Conclusions. 

II. DESCRIPTION OF NUMERICAL SIMULATIONS 

Before discussing how the simulations are carried out, we review the basic elements of the 
theory presented in Ref. ||^. We consider one atom at the origin, initially in the state s', in 
interaction with a gas of s atoms through a resonant ss' pp' dipole-dipole interaction V . 
We also allow for interaction among the unprimed atoms via an sp —>■ ps process mediated 
by another dipole-dipole potential U. This system is described by Eqs. (24a) and (24b) of 
Ref. p, which are 

N 

i do = Aao + J2 ^^Ck, (la) 

k=l 

N 

i Cfc = Vktto + J2 ^kici- (lb) 

1=1 

Here ao{t) is the amplitude of the state in which the atom at the origin is in state s' and 
all other atoms are in state s, while Ck{t) is the amplitude of the state in which the atom at 
the origin is in state p' and the atom at is in state p, while all the others remain in state 
s. The quantities Vk and Uki are the interaction potentials and A = + ep — e'^ — is the 
detuning from resonance. The potential Vk is of the dipole-dipole form 

\4 = -^P2(cos^,), (2) 

where 9k is the angle between and the magnetic field in the MOT used to confine the 
gas.Q Similarly Uki is of the form 



^The P2(cos0) angular dependence is also used, for instance, in Ref. p. In the absence of a 
magnetic field, each atom is in a multiplet of states and the angular dependences of each Vk and 
Uki are described by matrices. 
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[/H = -^P2(cOS^fc;), (3) 

for I 7^ k. We define Uu = 0, so that the restriction I ^ k is automatic in sums of the type 
appearing in Eq. (lb). Here /i is the dipole matrix element connecting the s and p states, 
and fi' is the dipole matrix element connecting the s' and p' states. The atoms are assumed 
to be sufficiently cold that during the time scale of interest they move only a very small 
fraction of their separation, and therefore and Uki can be taken to be independent of 
time. 

If we consider a column vector C (t) = (ao (t) , Ci (t) , C2 (t) , . . . , Cn (t))'^ , then clearly we 
can construct a matrix H such that the set of + 1 coupled differential equations (1) can be 
written in the form i C = HC . If we let {A^} and {6m} be the eigenvalues and eigenvectors 
(represented as column vectors) of the real Hermitian matrix if, respectively, and if we 
further define e„ to be the column vector of length + 1 with the nth entry set equal to one 
and all others set equal to zero, then we can write e„ = Z]m=o Pnmdm and 9m = Y.n=o Pmn^n, 
where Pnm = 6'^ ■ e„. The system starts off in the state described by Cq, so it follows that 

C{t) =exp{-iHt)eo 

N 

= PomGXp {-iXmt) 9m 

N N 

= PomPmn exp (-^A^^) e„. (4) 

m=0 71=0 

The quantity (t) is just the cq component of this expression, so 

N 

«o (t) = el ■C{t)= PomPmo exp {-iXmt) ■ (5) 

m=0 

Thus we see that, in this form, solving for ao at a given time becomes just a matter of 
diagonalizing the matrix H. To carry out the numerical simulations, we first generate a set of 
random positions for the unprimed atoms. For the results given in this paper these random 
positions are uniformly distributed, but it is just as easy to compute for a gas blob with a 
Gaussian density profile, as is the case in the experiments From these positions the 

potentials Vk and Uki and the matrix H are generated. We then diagonalize H and construct 
the quantity at a set of values of t using Eq. (|^). We do this for the same set of times 
for many random distributions of the unprimed atoms and thereby obtain an average of the 
quantities and |ao|^ at a given set of times over many configurations of the system. The 
signal measured by the experiments is then proportional to the average of 1 — |ao| . One 
can directly generate Fourier transforms from Eq. (|^), as we do in Section |V|. We can also 
obtain plots of the linewidth as a function of time by generating signal versus time curves 
for many values of the detuning. A, and then determining for each time value which of these 
signals is equal to one half the value of the resonance signal. 



^Results that directly compare calculations and experiments will be reported separately. 
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III. DESCRIPTION OF ANALYTICAL APPROXIMATION 



A. Justification 

In developing approximations, it is useful to cast the problem in the standard Green's 
function language [0, as discussed in Section 3 of Ref. |TT|. The Green's function of Eq. (lb) 



satisfies the equation 



N 

^Gin = 5/n + X! ^ImGmn, (6) 
m=l 

and then 

iO — L\ — 2^1^ ViLrimVm 

Actually, —iao is the 00 element of the Green's function for the entire set of + 1 equations. 
Analogously, then, the on-site Green's function for Eq. (lb) can be written in the form 

^ ~ J2lm UnlG im,[n]Umn ^ ^ 

where Gim,[n] is the Green's function for a medium in which the n atom is absent. It is 
important to keep in mind that these Green's function equations are exact. 

The Cayley approximation [|l5,|l^] keeps only the diagonal (/ = m) terms in the sum over 



/ and m in Eqs. (0) and (||). Further, it assumes that each Gii^[n] is independently distributed, 
and that there are sufficiently many sites that Gu^^n] and Gu have the same distribution. 

We choose to set uj = ia and work in Laplace space instead of Fourier space because 
this makes many of the integrals that follow more obviously convergent. Up to a point, 
the same analysis can be carried out with the Fourier analogs, but Laplace transforms are 
needed to solve the Cayley equations exactly in part B of this Section. We also work with 
the inverses of these Green's functions, setting ag (ia) = l//o (a) and Gkk (ict) = i/ fk («), 
so in the Cayley approximation we have 

N y2 

k=\ 
N tt2 

fk = a + Y.^- (9b) 

1=1 Ji 



B. U Process Only 

First we consider just the U process and solve Eq. (9b), which is to say we look at the 
sp exciton by itself. Note that fk is always real and positive. The distribution of fk is given 

by 
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\ \ 1=1, J I j I m^k,m=l 



(10) 



where for a uniform gas of volume Q, 

1 r ^ 



(11) 



Writing the deha function as a Fourier transform, we see that 

P{h) = r ^ I n P Um) dU e^^(^-") (exp (-^q ^] 

J -co m J ra^k,ra=l \ \ 1=1 Jl I 



dq 
2^ 



/ n PUm)df^—-l n d'ne^v[-iq^]. (12) 

mf^k,m=l " lytk,l=l \ Jl J 



Using the fact that (Hm Pm) (11/ -^/) = Dm PmEm, we can write 



P ifk) 



^ dq 
oo 27r 

oo 27r 
oo 27r 



oo 



ii fP (/-) rf/- ^ / ^r. exp (-iq^) 

k,m=l-' V ^"^ / 







MUk-a) 



dhP{f,)^Jd'nexp i-nj^ 



dfPif)^ I d'r {1-exp 



N-l 



-iq- 



(r) 
/ 



N-l 



(13) 



where in the last step we have simply changed variables from ri to r = ri — r^. Using the 
averaging methods detailed in Eqs. (8-10) and the Appendix of Ref. we see that 



P{f) 



where 



dq 
2^ 



Q{a) 



^iQif-o^) exp 



uJiqQ (a) 



dl_ 





P{f) 



and 



IGtt^/^ 2^ 
9^3 ^ 



(14) 



(15) 



(16) 



Using 



r 



dq „i 



iqA—y/ujB 



-oo 27r 

for A and real, we have 



1 B 

2v/^A3/2 



exp 



e(A) 



(17) 



6 



and so 



We also have 



or 



P{f) 



uQ (a) 



a 



,3/2 



exp 



"4(/-«). 



9{f-a). 



Q{a) 



uQ (a) 



,3/2 



exp 



M^Q^ (a) 
■4 (/-a) 



exp 



4a 



erfc 



(a) 
2v^ 



7. 



/ 



TT (a) 



3/2 
1 



exp 



u'Q^ {a) 



erfc 



(a) 
2v^ 



1 - (a) 



(18) 



(19) 



(20) 



(21) 



Thus we can solve numerically the transcendental relation for Q (a) given in Eq. ( p!9D and 
plug it into Eq. ( pl[ ) to find (1//)- It is important to note that Eq. (|T9|) defines Q as an 
analytic function of a in the half-plane where Re a > 0. Thus (1//) can be analytically 
continued, in the variable uj = ia, from the positive imaginary axis to the entire upper 
half-plane. 



C. U and V Processes in the Absence of A 



Now we allow the V process to be present, but we still ignore the detuning, A. We start 
with the equations 



N y2 

Jo = a + 2^ -Ti 

k=l 

N tt2 

fk = a+ 2^ — . 

i=i,i^k J I 

Repeating the analysis we did before, we have 

Af T/2 \ \ N 



which leads to the form 
Po (/o) = 



^ gig(/o-a) 
271 



dfP{f)^ I rfVexp 



—iq- 



V ir) 



f 



N 



(22a) 
(22b) 



(23) 



(24) 



At this point it is clear from our previous work that this equation leads to 



Poifo 



vQ (a) 



2A(/o-a) 



3/2 



exp 



v'Q^ (a) 



4(/o-a) 



(/o - «) 



(25) 



where 



and Q (a) is still given by Eq. (pj]), and depends on u. Therefore we obtain the result 

(1) r*P„(/.) 

\ /o / ^0 /o 

1 y/n vQ {a 



(26) 



a 2 a3/2 
which again can be analytically continued 



exp 



v'Q' (a) 
4a 



erfc 



vQ (a) 
2^^ 



(27) 



D. ?7 and V Processes in the Presence of A 

If we further consider the effect of a detuning, A, then we have the full set of Eqs. (9). 
The equation for fk is unchanged from the previous case, but because the equation for /o 
is now complex we have to consider a distribution in the two variables = Re/o and 
/g = Im/o. We have 

^0 (/o^ /oO = / (/(f - " - E f ) ^ (/o - A) ) n P ifi) dfi. (28) 

Clearly the delta function involving /q does not participate in the averaging or in the inte- 
grations over Therefore we can treat the part involving just as we did before, and the 
/q part will just be carried along. Hence we find 



Po (/(f) fo 



vQ (a) 



2v^ (/o« - a 



3/2 



exp 



v'Q' (a) 
4 (/o^ - a) 



(29) 



It follows, then, that 



vr vQ (a) 



a + iA 2 (a + iA) 



3/2 



exp 



v^Q^ {a) 
4 (a + iA) 



erfc 



vQ (a) 
27" + «A 



(30) 
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E. Other Approximations 



It is also possible to make approximations other than the one we have chosen. For 
example, instead of working with Eq. (9b), one can instead choose to work with 



fk 



a 



f ' 

1=1, J'' 



(31) 



which gives explicitly 



fk 



(32) 



The standard averaging procedure then leads to 



a 



u 



exp — erfc — 



u 



a 



Alternatively, we can replace 1/ fi by its average, obtaining 



(33) 



N 



fk 



a 



J. 



(34) 



which leads to 



a 



+ 



TT U 



U 



2 V^WfT'^^VocXf 



erfc 



u 



2a\\f 




1. 



(35) 



Eqs. (^Tj) and ( p^D correspond to a Hubbard band and a CPA approximation, respectively. 
We have examined these approximations and found that they give results that do not agree 
with our simulations as well as the results obtained by the Cayley approximation. This is 
not at all surprising, since in the Cayley equations (9) the on-site Green's functions for the 
U problem are allowed to vary from site to site, whereas in Eqs. (^) and ( P^ they are not. 



IV. RESULTS 

In this section we present and discuss the results of the numerical simulations described in 
Section |1| in comparison, when possible, with exact analytical results or with approximations 
that are good enough to be used instead of the simulations. We require that the graph of 
an approximate result should not entail much more numerical labor than the corresponding 
graph from the numerical simulations. Further, we omit comparisons of Laplace transforms, 
because two similar-looking Laplace transforms can correspond to distinctly different results 
in the time domain and to different spectral densities in the (real) frequency domain. We 
found that only the Cayley tree approximation for (oq), as developed in Section |T|, comes 



close to meeting all these criteria and can be considered as "good" for a broad range of 
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conditions, although in special limits simpler approximations suffice. We do not have, at 
this time, a good approximation for (|aoP). 

All the analytical results, exact and approximate, are for a single primed atom in an 
ideal, infinite gas of unprimed atoms with purely dipolar interactions. The simulations 
presented here were carried out for unprimed atoms randomly distributed in a cubic 
volume floating in vacuum, with a single primed atom at the center of the cube. One can 
change the position of the unprimed atom (to study surface effects), the shape of the volume, 
the boundary conditions, and the interaction. One can also consider a non-uniform gas blob, 
such as one with a Gaussian proffie, and put in some correlation between the positions of 
the atoms. Here we focus on the uniform, ideal g clS 8bS 8b standard. 

For compactness of presentation, all energy and inverse time variables are scaled by the 
effective strength of the ss' — pp' transition, v, defined in Eq. (0). We have then results as 
a function of time or frequency and of two parameters: A/v, where A is the detuning from 
resonance, and u/v, where u is the equivalent of v for the sp ps transition and is defined 
in Eq. (|T^). We display results for the following values of u\ 
u=0 

If the potential U is set equal to zero in Eq. (lb), then Eqs. (1) can be solved and 
averaged analytically as is done in Ref. [Q. This case therefore provides a convenient way to 
check to see if our numerical routines are working properly. It also provides a nice way to get 
a feeling for how many unprimed atoms need to be included and over how many realizations 
of the system we need to average to obtain good statistical results. In this case, of course, 
the sp exciton does not travel. 
u=v/4 

We report some results for this case in order to show the gradual transition from the 
non-propagating to the propagating exciton, and also because this value of u/v corresponds 
to the Rb case when there are, initially, a few s atoms in a gas of s' atoms (A^ <C A^'). In 
effect, one interchanges the roles of primed and unprimed atoms with respect to the case 
u = Av (see below), so that the s'p' exciton is now the one that propagates. 
u=v 

Resonant coupling to the excitons is nearly optimal in this case. In addition, for A = 
finding is equivalent to finding the on-site Green's function Goo for a gas of A^ + 1 atoms. 
Then — (l/7r)ImGoo (<^) = (l/7r)Reao iyj) is the spectral density for exciton propagation. 
This and other quantities related to exciton propagation can then be studied independently 
from the mechanism of resonant excitation. 
u=4v 

This value is chosen because it is close to the value u = 3.96v that pertains to the 
Rb experiments ll|,@,H. It also serves to illustrate, in general, what happens when exciton 
propagation is a large effect. 

A. Results for Npi /N' 

The four parts of Fig. |l] display the numerical simulations of 1 — {^|ao|^^, corresponding 
to the experimental signal Np/ /N', as a function of vt for u/v = 0,1/4,1 and 4, and enough 
values of the detuning A to give an idea of the resonance width in each case. The FWHM 
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at large t can be directly deduced from the intercepts of the graphs with the right edge of 
the figure, and is discussed more fully in subsection D. 

In addition, Fig. la shows the exact results for u = that were reported in Ref. [Q. 
Comparison with the numerical simulations, which were carried out for N = 100, shows 
that the infinite system is already well simulated, but that noticeable wiggles remain after 
averaging over 10000 realizations. The reason for this persistence of fluctuations is that 
most realizations do not look at all like the average. In fact, the variance with realization 
(not shown) is comparable to the average and in particular for small t the average is given 
by {y/n/2) vt, while each realization is quadratic in t, as already discussed in Ref. It can 
be seen from the graphs that the linear rise law, {^/2) vt, is independent both of A and 
of u and holds for a substantial range of vt. 



B. Results for (ag (t)) 



The four parts of Figs. |] and |^ display the numerical simulations of Re (oo) and Im (oq) , 
respectively, as a function of vt for w/f = 0, 1/4, 1, and 4. Recall that (ao) is the probability 
amplitude that resonant exciton creation has not occurred during the time t. The striking 
feature of these graphs is that, off- resonance, oscillations with an approximate period A 
persist for a long time, not only when the exciton does not propagate (m = 0), but even for 
u = 4f . 

The graph for u = can be obtained analytically by the same methods used in Ref. 
for computing 1 — (|ao|^). We give here only the following result for the Laplace transform: 



(ao (ia)) 



a + iA 2 (a + zA)^/'v^ 



exp 



4a (a + iA) 



erfc 



2Ja (a + iA) 



(36) 



From this, (ao (t)) has been obtained by expanding in powers of 1/a and inverting the 
Laplace transform term by term. As in Ref. |^ , various formulas involving special functions 
can also be obtained. For u = 0, the exact analytical results (not shown here) agree with 
the numerical simulations about as well as in Fig. la. The initial dependence on time is 
given by 



(37) 



TT 



(ao) ~ 1 vt — iAt 



Looking at the graphs, it seems that the initial decrease of Re (ao) depends on u, but this 
perception only indicates that the range of validity of the linear approximation (^) is small. 
On the other hand, the graphs of (||ao|^^ in Fig. || look linear and independent of u over 
a wide range of vt. Note also that |(ao)|^ ~ 1 — ^/nvt decreases initially twice as fast as 
<^|ao|^^ ~ 1 — {y/TT/2)vt. Of course, it must be true that the (mod) square of the average is 
smaller than or equal to the average of the square, but the above result shows that |(ao)|^ 
is not even a fair approximation to (|ao 
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C. Results for the Spectral Density 



Figs. § to 1^ present graphs of Re (ao {u>)) that correspond to the graphs of {a^ (t)) in 
Figs. ^ and ^ For u = v and A = 0, Re (oq (o;)) is the exciton spectral density (times 
vr); more generally, it is an excitation spectral density for the resonant process ss' —>■ pp'. 
Each Figure gives a comparison of the numerical simulation with the Cayley approximation, 
which in particular for u = reduces correctly to the exact result of Eq. (^6|), obtained from 
Eqs. (0) and ( jlQ]) by analytic continuation. We compare spectral densities, rather than 
time graphs, in part because obtaining (oq (t)) for the Cayley approximation is laborious, 
and in part because the comparison of spectral densities is instructive, as we now discuss. 

In the notation of Section |I|, the spectral density is obtained from Eq. (^) as 

Re ao {uj + ie) = Im ^ — /^om/^mo ^ ^^g^ 

in the limits e — > and N ^ oo. In practice it is difficult to go beyond N = 1000, and one 
must choose e to be fairly large to obtain a smooth plot for the spectral density where the 
eigenvalues are sparse, even after averaging over 10000 realizations. On the other hand, the 
distribution of eigenvalues for a random dipolar gas is highly structured near u = 0, and 
if e is too large this structure will be missed. Rather than spending much time to go to 
large and refine the numerical technique, we set N = 100, plot the results for e/v = 0.1 
and 0.01, and let the eye extrapolate to the correct limit. It is remarkable that the Cayley 
approximation agrees best with the smallest value of e/v, except for wiggles that are usually 
in regions of sparse eigenvalues. This is exactly how the exact spectral density is expected to 
compare with a calculation for finite A^. A finite e can also be viewed as a simulation of the 
decoherence introduced by atomic motions, or of losses to other channels, such as the decay 
of the excited atomic states due to blackbody radiation Thus, e is almost equivalent 
to the 7 of the Lorentzian broadening model of Ref. ||^, the difference being that in Ref. 
27 is added to u for the Fourier transform of Eq. (lb) only, while here ie is added to every 
UJ. As e increases, the highly structured spectral densities of the ideal frozen gas become 
increasingly similar to Lorentzian lines. Because of unavoidable losses and incoherences, the 
experimental spectral densities will resemble the results of the simulations for some finite e. 

For u = (Fig. 1) the spectral density is quite different from a Lorentzian. At A = 
there is a quasi-gap around uj = 0. As A increases the gap extends from u; = to = A. 

For u = v/4 (Fig. |^) the sharp features of the u = spectrum are rounded off and the 
gap is beginning to be filled, but is still recognizable. 

For u = V (Fig. ^ the spectrum at A = has the appearance of a single, asymmetric 
line. At A = f and A = 2f a shoulder on the left of the line is the remnant of the gap of 
Fig. |. 

For M = 4f (Fig. ^) there is a single line for all values of A, which becomes increasingly 
sharp as u increases. The lineshape is not Lorentzian, but even a small additional broadening 
will make it nearly so, as indicated by the calculation for e/v = 0.01. This sharp line leads 
to the damped sinusoidal oscillations seen in the u = Av graphs of Figs. |^ and |^. The general 
phenomenon of line narrowing can be discussed using the simple Lorentzian model for the 
exciton band, in which Eq. (lb) of Ref. [§| is replaced by 
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ick = Vkao- i'jCk. 



(39) 



Solving the equations in this limit, one finds a pole at 



which gives a decreasing width as u increases, since 7 is proportional to u. One cannot 
simply replace J2k by its average, which does not even exist. Instead, one can Fourier 
transform the average to find the time dependence 

exp (-iAt - , (41) 

where 

and V has been defined in Eq. (p6|). A similar argument was applied in Ref. [^] to find the 
time dependence of (|ao|^ 



From the comparison with the simulations, it is apparent that a shortcoming of the 
Cayley approximation is that it gives a symmetric spectral density at resonance (A = 0), 
while the correct spectral density is slightly asymmetric, except for m = 0. Correspondingly, 
(ao it)) is not purely real, as can be seen in Fig. |^. For A 7^ 0, the spectral density at a; < A 
is underestimated by the Cayley approximation. 

Clearly, the asymmetry of the spectral density at A = comes from the terms with 
/ 7^ m in the denominator of Eq. (|^), which are neglected in the Cayley approximation. It is 
surprising that these terms, which correspond to exciton propagation other than back and 
forth on the Cayley tree, have such small effect on (ao) at A = 0. At finite A, the underes- 
timate of Re (flo iyj)) for < A is probably due to the neglect of these same terms. It seems 
complicated to include these terms, even approximately, in the Cayley tree formulation, but 
we have already shown how this can be done in the "CPA" approximation [O . 



D. Results for the Resonance Linewidth 

The time-dependent resonance lineshape is obtained from the ordinates of the curves in 
Fig. H at a given vt. In general the lineshape is not Lorentzian, but we never find a split 
line with a minimum at resonance as reported in Ref. . Thus it is fair to characterize the 
line by its FWHM. This cannot be directly compared with an experimental width, which is 
usually the FWHM of a Lorentzian fitted to a noisy experimental line, but it is indicative 
of trends. For a more accurate comparison with experiment, one should of course convolute 
the calculated lineshape with the broadening from other effects. 

Fig. 1^ shows one half of the FWHM, w, as a function of vt for u = 4v. One sees the 
typical decrease of w from a universal behavior at small vt towards an asymptotic value, 
which is almost reached at vt = 15. The t~^ law, which results simply from the uncertainty 
principle, has already been discussed in Ref. ||^, where on the basis of a small A expansion 
it was estimated that w = \/l2/t for small t. 
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The near-asymptotic value of w at vt = 15 is shown in Fig. ^ as a function of u/v. 
Overall, the graph is roughly linear except at small u/v, although a careful examination 
shows that there are really two linear regimes with somewhat different slopes that meet at 
u/v ^ 3.5. There is an unexpected minimum around u = v/A, which is not a numerical 
artifact, but can be seen directly from Fig. |1| by comparing the plot for u = v/A with those 
for u = and u = v. What happens is that the resonance line always grows taller and wider 
as u increases, but for small u its height (the signal at A = 0) grows faster than its width. 



V. CONCLUSIONS 

In this paper, we report results of ongoing work on the general question: what are the 
electrodynamics and the optics of a gas of excited atoms interacting via dipolar fluctuations, 
for times short compared to the atom-atom collision time? We are considering in particular 
the time evolution of a resonant ss' — > pp' process when N atoms in state s and N' in state s' 
are initially present, and all atoms are randomly distributed. At least one of these two initial 
states must be an excited state for the process to occur resonantly. Other resonant processes 



can be mapped on this prototype. For instance, the process pp — > ss' discussed in Refs. fTS 
and 1^ is handled by setting s' = s and reversing the roles of s and p. If the atoms are in 
Rydberg states then they have known energies and dipole moment matrix elements and the 
dipole-dipole interaction is dominant, so there are no adjustable parameters in the theory. 
More importantly for the experiments, with a frozen gas of Rydberg atoms one can easily 
do time-resolved spectroscopy because the evolution time scale is measured in microseconds. 
However, for the theory it does not matter that the atoms are of Rb or Cs, since other atoms 
in excited states would behave similarly on the corresponding time scale. In the language of 
the solid-state and chemical physics literatures, our problem can be described as resonant 
energy transfer to Frenkel excitons. 

We have carried out detailed and extensive calculations, truly "computer experiments" , 
for a resonant ss' pp' process when a single s' atom is initially present in a frozen gas of 
s atoms. Neglecting the propagation of the s'p' exciton to keep things simple, we consider 
in effect the resonant excitation and propagation of an sp exciton. The results apply, in 
practice, whenever only a few excitons of a given type are created, so that exciton-exciton 
interactions are negligible. 

We can compute anything we wish about this system, but in this paper we report rep- 
resentative results for two quantities only: the average probability amplitude that no exci- 
ton has been created, (ao), and the average probability that an exciton has been created, 
1 — ^|ao|^^. One striking result is that the initial transition rate for this process is indepen- 
dent of exciton propagation and of detuning from resonance for surprisingly long times, and 
can thus be explicitly evaluated. Another notable feature is that the transition amplitude 
at detuning A shows persistent oscillations with period A, in particular when exciton prop- 
agation is a strong effect, as for u = 4v. Correspondingly, the spectral density for exciton 
creation is strongly peaked at A, and becomes narrower as u increases. Thus, the uni- 
versal phenomenon of line narrowing actually simplifies the problem in the experimentally 
interesting limit of strong coupling (m t>). 
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With a wealth of numerical results at our disposal, we have attempted to develop approx- 
imate formulas that are not simply fits to a particular numerical result. We have obtained 
good approximations for (oq), but not yet for (jaof'^, which is notoriously more difficult 

to obtain than (ao) and quite different from |(ao)|^. To develop approximation schemes we 
have drawn, and will continue to draw, from the literature on wave functions in random 
systems, and in particular on earlier work on excitons in disordered solids and liquids |]T3 . 
The best results, in our case, are given by the Cayley tree approximation. 

Most of the work on random systems deals with on-site disorder or short-range inter- 
actions. On the contrary, the dipolar interaction U (r) ~ is long range and has the 
distribution P (U) ~ for a random distribution of r. The distribution of the cumulative 
variable U, where = J2k U"^ {'^k)i plays a major role in the theory. It is given by 



^(W) = ^77^exp --- , (43) 




where u is given by the appropriate equivalent of Eq. Note that P (JA) is not Gaussian 

and the central limit theorem does not apply, because P (U) does not have an r.m.s. value and 
is not even normalizable. Therefore it is not surprising that, as shown by the simulations, 
the distribution of eigenvalues and other statistical properties of the matrix Uki are very 
different from those of a standard random matrix. This implies that random system results 
obtainable by random matrix theory are not applicable to dipolar interactions in an ideal 
gas. Luckily, the distribution in Eq. ( ^31 ) makes the Cayley tree equations for (ao) exactly 
soluble in this case. 
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^The angular dependence of U modifies the value of u, leaving unchanged the functional form of 
Eq. (|3|). See Ref. §, Appendix A. 
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FIG. 1. Plots of 1 - (jao {t)f) loTu = 0,u = 0.25v, u = V, and u = 4:V for several values of 
the detuning A, as shown on each graph. For u = 0, the exact result is shown for comparison. 
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FIG. 2. Plots of Re (ao [t)) ioi u = 0, u = 0.25v, u = v, and u = 4v for three values of the 
detuning A. 
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FIG. 3. Plots of Im (ao (t)), for the same values of u and A as in Fig. |2|. 
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FIG. 4. Plots of Re (oq (uj)) for n = and A = 0, A = v, and A = 2v. Exact results are 
compared with numerical simulations of Re (ao (a; + ie)) for two values of e. 
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FIG. 5. Plots of Re (ao {to)) for u = 0.25v and A = 0, A = v, and A = 2v. 



40 




42 



6 



Cayley Result 
£/v=0.01 

£/V=0.1 



o 
(0 

CD 
DC 



D) 

CO 
i_ 

CD 

> 

< 







-5 -4 -3 -2 



■1 1 
co/v 



u=v 
A=2v 





3 4 5 



43 



FIG. 6. As in Fig. |5|, but for u = v. 
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FIG. 7. As in Fig. |, but for 
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FIG. 8. A plot of the width as a function of time for 
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FIG. 9. A plot of the width at vt = 15 as a function of the ratio u/v. 
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